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We analyze the dynamics of dissipative solitons in silicon on insulator waveguides embedded in a gain 
medium. The optical propagation is modeled through a cubic Ginzburg-Landau equation for the field 
envelope coupled with an ordinary differential equation accounting for the generation of free carriers owing to 
two-photon absorption. Our numerical simulations clearly indicate that dissipative solitons accelerate due to 
the carrier-induced index change and experience a considerable blue-shift, which is mainly hampered by the 
gain dispersion of the active material. Numerical results are fully explained by analytical predictions based on 
soliton perturbation theory. © 2013 Optical Society of America 
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Silicon photonics has attracted considerable attention 
among researchers owing to its potential applications 
ranging from optical interconnection to bio-sensing. In 
the last few years, silicon-on-insulator (SOI) technology 
has rapidly developed into a well-established photonic 
platform [1] . The tight confinement of the optical mode 
and the inherently large bulk nonlinearity of Si tremen- 
dously enhance the nonlinear dynamics [2]. For near- 
infrared wavelengths in the range 1 fxm < Xq < 2.2 /im, 
two photon absorption (2PA) is the leading loss mecha- 
nism and limits the spectral broadening due to self-phase 
modulation (SPM) [3,4]. As a consequence of 2PA, elec- 
trons are excited to the conduction band, absorb light 
and affect the pulse dynamics by modifying the refrac- 
tive index of silicon [5] . Loss mechanisms are restricted 
in silicon-organic hybrid slot waveguides, which can be 
exploited for all-optical high-speed signal processing [6] . 
Alternatively, loss can be overcome by embedding active 
materials in the design of SOI devices. Recently, some 
amplification schemes based on III-V semiconductors, 
rare-earth-ion-doped dielectric thin films and erbium- 
doped waveguides have been proposed and practically 
realized [7-9]. 

In this Letter, we describe for the first time to our 
knowledge the nonlinear dynamics of dissipative solitons 
(DSs) in amplifying SOI devices. DSs are stationary lo- 
calized structures of open nonlinear systems far from 
equilibrium that can be observed in several contexts [10]. 
Although the model considered in our analysis is gen- 
eral and can be applied to any silicon-based amplifying 
setup, we specialize our calculations to a SOI waveg- 
uide embedded in Er-dopcd amorphous aluminium ox- 
ide (Al203:Er + ). Other gain schemes involving the use 
of semiconductor active materials can also be consid- 
ered and gain dispersion can be reduced accordingly. For 
the representative structure displayed in Fig. 1, the gain 
bandwidth is of the order of 100 nm around the carrier 
wavelength A ~ 1540 nm. For this waveguide, the second 
order GVD coefficient and the effective area at Ao = 1550 



nm are calculated to be /3 2 — —2 ps 2 /m, A c g ~ 0.145 
fim 2 , respectively. The proposed SOI waveguide is fabri- 
cated along the [110] direction and on the [110] x [001] 
surface, so that quasi-TM modes do not experience stim- 
ulated Raman scattering (SRS) [5]. The propagation of 
an optical pulse with envelope u(z, t) and carrier fre- 
quency loq in the proposed photonic structure is governed 
by a complex Ginzburg-Landau (GL) equation, 

id^u — isgn(/3 2 )<9 2 u + iau + (1 + iK)\u\ 2 u 

+(i/2- n)(f> c u - i(g + : g 2 d 2 )u = 0, (1) 
coupled with an ordinary differential equation account- 
ing for the 2PA- induced free-carrier dynamics d<j) c / dr = 



e\u\ 



Eq. (1) is written in dimcnsionlcss units, 



where the time (t) and the longitudinal spatial variable 
(z) are normalized to the initial pulse width (t = tor) 
and to the dispersion length (z — £Ld with Ld = 
tQ/\(3 2 {^o)\): respectively. The envelope amplitude (A) 
is rescaled to A = u^/Po, where Pq = |/32(wo)|/(io7fl)j 
7ii = k n 2 /A cS with n 2 ~ (4± 1.5) x 10~ 18 m 2 /W is the 
Kerr nonlinear coefficient of bulk silicon. The bulk 2PA 
coefficient is /?tpa — 8 x 10~ 12 m/W and its correspond- 
ing effective waveguide counterpart 7/ = /?TPA/(2Arff ) 
is rescaled to the Kerr coefficient so that K = ji/-fji — 
/3tpaAo / '(47rn2) . The linear loss coefficient (a;) is renor- 
malized to the dispersion length (a = c^Lrj) and can 
be neglected for short propagation in the linear trans- 
parency spectral window of silicon 1 /xm< Ao <10 /im, 
where 2PA dominates (if Ao < 2.2 ^m). The density of 
free-carriers (FCs) N c generated through 2PA is normal- 
ized so that 0c = <tN c Ld, where a ~ 1.45 x 10 _21 m 2 . 
FCs are responsible for free-carrier dispersion (FCD) reg- 
ulated by the parameter 6 = /?TPA|/32|o'/(2/iwo^4cff i o7^) 
[11] and free-carrier absorption (FCA) depending on the 
parameter /i = 2irk c /(a\o), where k c ~ 1.35 x 10~ 27 m 3 
[12]. t c is the characteristic FC recombination time (of 
the order of ns) and normalized as t c = to/t c , which 
we neglect in our calculations since we focus our analy- 
sis on ultrashort pulses with time duration of the order 
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Fig. 1. Sketch of a SOI waveguide with lateral dimensions 
h — w — 525 nm surrounded by Al203:Er + and its disper- 
sion property. The solid red line represents the group velocity 
dispersion (GVD) of the quasi-TM mode, whose spatial pro- 
file at Ao = 1550 nm is depicted in the inset below. 



to ~ 100 fs. The amplifying medium is characterized by a 
gain coefficient G, which in our dimensionless equations 
is rescaled to the dispersion length (g = GLd), and a 
dephasing time T2, which is related to the dimensionless 
gain dispersion coefficient through cj2 = g(T2/to) 2 . 

The presence of gain and 2PA in the photonic struc- 
ture considered in our calculations implies that, in gen- 
eral, energy is not conserved. The GL equation does not 
support solitons in a strict mathematical sense since it 
is not integrable by the inverse scattering transform. In 
contrast to Kerr solitons in conservative systems, which 
form continuous families of localized solutions, DSs are 
formed under a dynamical equilibrium involving non- 
trivial internal energy flows. DSs are associated with cer- 
tain discrete parameters of the GL equation that satisfy 
the energy balance condition. In absence of linear loss 
(i.e. a = 0) and free-carriers (i.e. cf> c = 0), DSs can be 
found as u(£,t) = Uo[sech(r)T)] 1 ~ z P e lT ^ , where the exact 
values of uo,r),/3,T are fixed by the physical parame- 
ters of the system g,g2 and K [13]. This soliton suffers 
from inherent core and background instabilities and can 
be stabilized by introducing higher order nonlinearitics 
or by coupling the system to a passive waveguide [14]. 
However, this task goes beyond the scope of the present 
work, where we aim at understanding the effect of FCs on 
the DS dynamics. In order to grasp the fundamental ef- 
fects induced by FC dynamics, we develop a soliton per- 
turbativc analysis [15], approximating the GL equation 
as a perturbed nonlinear Schrodinger equation (NLSE): 



ie(u), where we explicitly con- 



sider the case of anomalous dispersion (/?2 < 0) and e(u) 
includes the coupling to FCs, 2PA, linear gain and its 
dispersion: e = (g — K\u\ 2 — 4> c /2 — i^<p c + g2d 2 )u. The 
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Fig. 2. (a) Asymmetric spectral broadening of the NLSE soli- 
ton (uin — sechr) with time duration to = 40 fs for g = 1 and 
<?2 = 0. (a) DS amplitude tin as a function of linear gain g for 
several dephasing times (t2 = T2/io = 0, 0.2, 0.5) at fixed K. 
(c) Carrier-induced blue-shift of DS (utn = tto[sec/i(r;r)] 1_1 ' 3 ) 
for the same parameters of (a), where uq = ^3g/(2K), 
V = \f I FfP and /3 = 3/(2#) ~ a/[3/(2^)] 2 + 2. Solid red 
and blue dashed lines indicate the output and input power 
spectra, respectively, (d) ^-dependent DS frequency shift 
(AS = Q/(2n)) for (72 = 0. Perturbative predictions with 
constant (solid red line) and z-dependent (green dotted line) 
peak amplitude un are shown. The solid red dots indicate 
numerical findings. 



perturbative theory is developed by making the Ansatz: 
<Cr) = UoiOisechMOlT-TpiO]}} 1 ^ * 

where the parameters Uo, rj, t p , <j), fl are now assumed to 
depend on £ and e(u) is considered as a small perturba- 
tion depending on u, u* and their derivatives. The evo- 
lution dynamics of the soliton parameters over distance 
can then be predicted using the variational method [13], 
which leads to a set of coupled differential equations for 
the soliton parameters. The evolution of pulse energy 
(E), frequency (f2) and temporal (r p ) shifts are given by 

^ = -2E^-eulE + 9 2^y (3) 



dn 

drp 



15 



( f i + P/2)6ut--g 2 (l + P 2 W 



(l-2g 2 m 



72 



-9E 



(4) 
(5) 



The perturbative analysis reveals that FCs induce fre- 
quency blue-shift and temporal acceleration of DSs, both 
effects being hampered by the gain dispersion that lim- 
its blue-shifting within the amplifying frequency win- 
dow of the active material. For small gain dispersion, 
Eq. (4) is solved with the initial condition 17(0) = 0, 
obtaining = /£, where / = 8([i + /3/2)^/15. 

This expression predicts that DSs experience a spec- 
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Fig. 3. (a,b) Input (dashed blue line) and output (£ = 5, full 
red line) power spectrum of a dissipative soliton with carrier 
wavelength Ao = 1550 nm and time duration to = 200 fs for 
g = 1,<72 = 0.04 as a function of (a) dimensionless frequency 
and (b) physical wavelength, (c) Input (dashed blue line) and 
output (£ = 5, full red line) temporal DS profile, (d) Intensity 
counterplot of the spatio-temporal evolution of an accelerated 
dissipative soliton. 



tral blue-shift proportional to the propagation distance 
through a rate (/) that basically depends on the FC 
density. The equation of the temporal shift turns out 
to be 7p(£) = -(af + /C 2 /2), where a = 70E 2 /72. The 
expression suggests during the propagation DS is accel- 
erated under the influence of FCs analogously to the 
recently studied case of gas-filled hollow-core photonic 
crystal fibers [16]. 

In Fig 2a we show the conventional asymmetric spec- 
tral broadening of NLSE solitons [17]. In Fig. 2b we plot 
the DS amplitude (uq) as a function of linear gain for 
several dephasing times (t 2 = T2A0 = 0,0.2,0.5). As 
shown in Fig. 2c, even though DSs of the cubic GL 
equation arc unstable, they experience a considerable 
carrier-induced blue-shift maintaining their shape over a 
propagation distance of the order of millimeters. Indeed, 
for DSs, 2PA is exactly compensated by the linear gain, 
while NLSE solitons experience an imbalanced dynam- 
ical evolution due to 2PA and amplification. In Fig. 2d 
we compare the perturbatively predicted FC-induced fre- 
quency blue-shift (full blue and dashed green lines) with 
numerical results (solid red dots), finding that for the ini- 
tial part of pulse propagation, perturbative predictions 
nicely match with the numerical findings. The dimen- 
sionless frequency shift is calculated to be Ad ~ 0.16 
at £ = 2, corresponding in physical units to A A ~ 35 
nm. The solid red line in Fig. 2d represents the pre- 
dicted frequency shift calculated through the perturba- 
tion analysis by approximating the DS amplitude to re- 
main constant, whereas the green dotted line represents 
the perturbative prediction considering ^-dependent DS 



amplitudes. From the very beginning of the propaga- 
tion, the pulse shape is perturbed by FCD and FCA 
and hence the approximate analytic treatment fails for 
long propagation distances. Fig. 3 displays the effect of 
gain dispersion over the pulse propagation. Due to the fi- 
nite bandwidth of the amplifying medium, the frequency 
shift is hampered after some saturation frequency (f2 sa t), 
as clearly predicted by our perturbative analysis: f2(£) — 
fiaatCl-e-*), where n sat = {ii + fi/2)0E 2 /[lQg 2 (l + fi 2 )\ 
and p = (4/3)52(1 + /? 2 )?7 2 - In Figs. 3a,b we depict the 
input and output spectral power of a DS in presence of 
gain dispersion indicating that blue-shift is reduced, as 
predicted by the theory. In Fig. 3c we show the temporal 
profile of a propagating pulse at £ = 5 where the inherent 
background instability of DSs starts to affect the pulse 
propagation. The spatio-temporal pulse evolution is dis- 
played in Fig. 3d, where the pulse acceleration is reduced 
by the effect of gain dispersion. 

In conclusion, we have provided a complete theoretical 
analysis of the the carrier-induced DS dynamics in Si- 
based waveguides embedded in an amplifying medium. 
FCs generated through 2PA affect the refractive index 
of the medium and lead to a considerable self-frequency 
blue-shift. We have derived analytical predictions for the 
self-frequency blue-shift and temporal evolution based 
on soliton perturbative theory. We also examined the 
fully realistic condition where gain dispersion hampers 
the continuous spectral blue-shifting, which is still ob- 
served. 
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